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INTRODUCTION 

Problem  Statement  and  Scope 

The  objective  of  this  thesis  is  to  present  an  ultimate  strength 
analysis  of  composite  beams  with  web  openings.  A  composite  beam  is 
defined  as  a  steel  W  shape  acting  together  with  a  concrete  slab  to 
resist  transverse  loads.  An  opening  located  in  the  web  of  the  steel 
section  is  usually  introduced  to  permit  the  passage  of  utility  ducts 
and  piping.   Figures  1  and  2  show  elevation  and  cross  section  views  of 
a  composite  beam  with  a  web  opening. 

The  analysis  is  limited  in  scope  by  the  physical  characteristics 
of  the  beam,  and  the  type  of  failure  assumed  at  the  opening.   The  slab 
thickness  is  limited  to  the  range  of  values  normally  encountered  in 
practice,  and  the  slab  width  is  taken  to  be  the  effective  width,  which 
is  determined  in  the  usual  manner  (11).  A  sufficient  number  of  shear 
connectors  are  assumed  to  be  present  so  that  full  composite  action  is 
attained.  The  opening  is  limited  to  a  rectangular  shape,  which  can  be 
located  anywhere  on  the  span,  and  can  be  concentric  (mid-depth  of  opening 
coincides  with  mid-depth  of  steel  shape)  or  eccentric.   Only  unreinforced 
openings  are  considered.   Failure  is  limited  to  yielding  only,  i.e., 
buckling  and  instability  failures  are  not  considered. 

Review  of  Previous  Ultimate  Strength  Analyses 

In  the  past  decade  a  number  of  investigators  have  developed  ultimate 
strength  analyses  of  non-composite  beams  with  rectangular  web  openings. 
All  of  these  analyses  lead  to  the  development  of  an  interaction  diagram 
which  shows  the  relationship  between  moment  and  shear  acting  at  an 
opening  at  failure.   Several  basic  assumptions  are  common  to  these 


analyses.  A  failure  mechanism  is  assumed  to  form  with  plastic  hinges 
located  at  the  sections  above  and  below  each  edge  of  the  opening. 
Failure  due  to  instability  is  not  considered.   Equilibrium  conditions 
are  satisfied.  Yielding  occurs  in  the  flanges  due  to  tension  or 
compression,  and  yielding  in  the  web  due  to  combined  shear  and  normal 
stresses  follows  von  Mises  yield  criterion  (10) .  The  presence  of  shear 
causes  secondary  moments  in  the  top  and  bottom  sections.  None  of  the 
analyses  take  into  consideration  the  beneficial  effect  of  strain 
hardening. 

The  first  analysis,  which  was  concerned  with  concentric  openings 
with  no  reinforcement,  was  developed  by  Bower  (1).   The  possibility  of 
the  web  and  flanges  having  different  yield  stresses  was  provided  for 
in  this  analysis.   The  shear  force  was  applied  only  to  that  portion  of 
the  web  which  was  also  assigned  the  secondary  moment.  Later,  in  dealing 
with  the  same  case,  Redwood  chose  to  have  the  same  yield  stress  through- 
out the  section,  and  also  assigned  the  shear  force  uniformly  along  the 
total  depth  of  the  remaining  web  (7).  Redwood's  revisions  were  incor- 
porated into  subsequent  analyses  of  concentric  reinforced  openings  by 
Congdon  and  Redwood  (2),  eccentric  unreinforced  openings  by  both  Frost  (4) 
and  Richard  (8),  and  the  most  general  case  of  eccentric  reinforced  open- 
ings by  Wang  (12). 

New  insight  for  the  analysis  of  beams  with  web  openings  was  presented 
in  a  report  by  McCormick  (6).  By  the  use  of  two  new  concepts,  McCormick 
developed  a  much  simpler  analysis  than  any  of  those  previously  presented. 
One  of  these  concepts  is  to  assign  a  moment  due  to  eccentricity,  M  ,  in 
the  larger  tee  section  to  represent  the  stresses  in  that  section.  As  in 
previous  analyses,  the  shear  force  was  assigned  to  the  full  web  stub 


length,  but  in  applying  von  Mises  criterion  the  web  thickness  was 
reduced  according  to  the  value  of  shear  present,  so  that  the  effect 
of  the  shear  stress  can  be  ignored  throughout  the  remainder  of  the 
calculations.   Because  of  these  new  concepts — introduction  of  M  and 
reduction  of  the  web  thickness  for  shear — axial  forces  and  moments, 
instead  of  stress  blocks,  were  used  in  a  statical  method  for  a  lower 
bound  approach  which  leads  to  a  simpler  analysis. 

A  comparison  between  Redwood's  and  McCormick's  analyses  was  made 
by  Scritchfield,  who  concluded  that  "McCormick's  method  of  analysis 
was  found  to  be  better  suited  for  extension  to  the  eccentric  case"  (9) . 
Scritchfield  applied  McCormick's  method  to  the  case  of  eccentric  un- 
reinforced  web  openings  by  the  use  of  a  computer  program,  which  when 
compared  with  earlier  programs  using  Redwood's  method,  gave  the  same 
results.   It  was  also  proved  that  the  points  of  contraf lexure  are  at 
the  center  of  the  opening. 

The  only  material  reviewed  pertaining  to  ultimate  strength  analysis 
of  composite  beams  with  web  openings  was  that  found  in  McCormick's 
report  (6).   In  the  report,  McCormick  performs  an  analysis  of  a  specific 
composite  beam  with  known  dimensions  and  material  properties,  having 
two  circular  web  openings  with  varying  types  of  reinforcement.   The- 
assignment  of  internal  forces  is  carried  out  in  a  manner  similar  to  that 
used  for  non-composite  beams.   The  concrete  slab  is  assumed  to  carry 
no  shear.  An  equivalent  rectangular  opening  having  a  depth  of  0.9D  and 
a  width  of  0.A5D,  where  D  is  the  diameter  of  the  circular  opening,  is 
assumed  for  the  failure  mode  consisting  of  a  four  hinge  mechanism  at 
one  opening.   McCormick  also  assumes  a  constant  distance  between  the 
axial  forces  in  the  top  and  bottom  tees  instead  of  determining  this 
distance  from  beam  properties  for  each  value  of  total  shear  force. 


The  analysis  presented  in  this  thesis  has  many  assumptions  in 
common  with  McCormick's  analysis,  but  is  developed  for  general  beam 
geometry  and  material  properties,  and  for  a  single  rectangular  opening 
of  any  practical  depth,  width,  and  position. 


ULTIMATE  STRENGTH  ANALYSIS 
Assumptions 

The  ultimate  strength  analysis  is  based  on  the  following  assumptions: 

1.  The  compressive  strength  of  the  concrete  in  bending  is  assumed  to 
be  0.85  f^  and  the  Whitney  stress  block  is  used. 

2.  The  tensile  strength  of  the  concrete  is  neglected;  therefore 
yielding  in  the  concrete  is  by  compression  only. 

3.  Yielding  in  the  steel  flanges  is  by  compression  or  tension  only. 

4.  Shear,  which  causes  secondary  bending  in  the  sections  above  and 
below  the  opening,  is  carried  in  the  web  only,  and  is  uniformly 
distributed. 

5.  Yielding  in  the  web  of  the  steel  section  due  to  combined  shear 
and  normal  stresses  follows  von  Mises  yield  criterion. 

6.  Equilibrium  is  satisfied. 

7.  Points  of  contraflexure  occur  at  the  midpoints  of  the  sections 
above  and  below  the  opening. 

8.  Failure  occurs  by  the  formation  of  a  mechanism  with  hinges 

at  sections  above  and  below  the  edges  of  the  opening.   (Fig.  3). 

9.  The  possibility  of  failure  due  to  instability  and  the  beneficial 
effects  of  strain  hardening  are  not  considered. 

Outline  of  Solution 

The  solution  is  divided  into  two  parts,  designated  Case  I  and  Case 

II.   Case  I  is  called  the  low  shear  case,  during  which  all  of  the  total 

shear  force,  V,  assigned  to  the  beam  is  carried  by  the  top  tee,  i.e.,  the 

shear  in  the  top  tee,  VT,  equals  the  total  shear  V.   Because  no  shear 

force  is  assigned  to  the  bottom  tee  in  Case  I,  the  capacity  of  the  bottom 

tee  is  used  solely  for  the  axial  force  Pfi,  which,  when  combined  with  an 

equal  force  in  the  slab,  gives  the  primary  moment,  P_d  . 

B  c 

A  special  situation  to  consider  at  the  outset  of  Case  I  is  that  of 
pure  bending,  i.e.  V  =  0,  (Fig.  4a).  The  total  capacity  of  the  top  tee 
is  assigned  to  the  axial  force  PT,  which,  when  combined  with  an  equal 


force  in  the  slab,  results  in  the  moment  due  to  eccentricity,  M  =  P  d  . 

The  moment  capacity  at  the  centerline  of  the  opening  is  the  sum  of  the 

primary  moment,  PRd  ,  and  M  . 

When  the  shear  force  in  Case  I  is  non-zero,  the  web  thickness,  t  , 

w' 

of  the  top  tee  is  reduced  to  w  according  to  von  Mises  yield  criterion, 

so  that  all  the  fibers  in  the  reduced  steel  section  will  be  at  the  yield 

stress.  A  secondary  moment  due  to  shear,  M^  =  V  a  is  induced  in  the  top 

tee  (Fig.  4b).   This  causes  a  reduction  in  P„,  and  likewise  in  M  .  The 

l  e 

total  moment  capacity  at  the  centerline  of  the  opening  is  still  the  sum 

of  the  primary  moment,  P^d  ,  and  M  .   The  upper  limit  of  Case  I  is  reached 

dc       e 

when  the  total  top  tee  is  yielded  due  to  V  and  M^  ,  so  that  M  is  equal 
to  zero. 

Case  II  (Fig.  4c)  is  called  the  high  shear  case  during  which  part  of 
the  total  shear  goes  to  the  top  tee  and  the  rest  goes  to  the  bottom  tee. 
The  amount  of  the  total  shear  assigned  to  the  top  tee  is  governed  by  the 
capacity  of  the  top  tee  section  for  V  and  M^  =  V  a.   The  amount  of  shear 
remaining  when  this  capacity  is  reached  is  the  shear  assigned  to  the  bottom 
tee,  Vg.  With  shear  present,  the  web  thickness  of  the  bottom  tee  is  reduced 
to  w  ,  and  a  secondary  moment  due  to  shear,  My  =  VRa,  is  induced.   The 
axial  force  P  is  assigned  to  that  portion  of  the  bottom  tee  not  used  for 

VB  or  ^Vfi*   Ttie  force  Pb'  alonS  wittl  an  equal  force  in  the  concrete  slab, 
gives  the  primary  moment,  which  is  the  total  moment  capacity  at  the  center- 
line  of  the  opening,  because  M  is  zero  throughout  Case  II. 

Development  of  Basic  Equations 

Reference  Values.  At  the  outset,  a  number  of  reference  values  are 
defined.   The  length  of  the  web  stubs  above  and  below  the  opening  are 
(Fig.  2) 


sT  =  %  d-e-h-t  (1) 

sB  =  %  d+e-h-t  (2) 

The  shear  capacities  of  the  top  and  bottom  web  stubs  by  definition  are 
STtwFy 

V       ^-  ^> 

s-,t  F 
v   «  _J_5L^  (4) 

yB       ^  (4) 

From  Fig.  5a,  the  shear  capacity  of  the  web  without  the  opening  (the 
gross  web  area)  is 

(d-2t)t  F 

The  total  plastic  moment  of  the  gross  composite  section,  Mp  ,  is  the 

final  reference  value  required.   Two  expressions  for  >L,  are  possible 

depending  on  the  location  of  the  plastic  neutral  axis,  NAp  of  the  gross 

composite  section.   To  determine  where  this  neutral  axis  is,  a  comparison 

is  made  between  the  total  axial  force  capacity  of  the  concrete  slab 

P   =  b  CF  (6) 

yc    c  c  KOJ 

and  the  total  axial  force  capacity  of  the  gross  steel  section 

Pys  -  (tw(d-2t)  +  2bt)Fy  (7) 

If  P   is  greater  than  P  ,  then  the  NAp  is  in  the  concrete  slab  as 

shown  in  Fig.  5a.   The  thickness  of  concrete  used  to  give  a  force  in  the 

concrete  slab  equal  to  that  of  the  steel  section  is  given  by 

P 
ys 

cPs  =  bT  (8) 


c  c 


This  is  the  thickness  of  the  concrete  above  the  NA^;  the  concrete  below 
the  NAp  is  disregarded  or  "thrown  away"  because  it  is  in  tension.   The 


value  of  the  total  plastic  moment  is  found  by  summing  the  moments  about 
the  NAp  resulting  in 

*?c   "  ^VpsX  +  C%d  +  c  -  cps)Pys  (9) 

If  P   is  less  than  P  ,  the  NA,,  is  in  the  top  steel  flange  as  in 
yc  ys       T 

Fig.  5b.   To  find  its  location,  a  thickness  t  is  assigned  to  the  portion 
of  the  flange  which  is  in  tension  below  the  NAp.   By  setting  the  forces 
above  and  below  the  NA_  equal  to  each  other,  the  value  of  t  is 

b  cF  -  t  (d-2t)F 

y 

Now  by  summing  moments  about  the  NAp,  the  total  plastic  moment  is 
Mpc  =  bcc(%c  +  t  -  tt)Fc  +  [tw(d-2t)(%d  -  t  +  tt) 

+  %b(t-tt)2  +  %btfc2  +  bt(d  -  ■—■  +  tt)]F  (11) 

Low  Shear  Solution.   The  following  discussion  of  the  analysis  is 
divided  into  two  major  parts:   Case  I  being  the  low  shear  case  and  Case  II 
being  the  high  shear  case.   In  Case  I,  the  total  shear  force  is  applied 
to  the  top  tee,  i.e.  V  =  V.   In  assigning  this  shear  force  to  the  web, 
a  portion  of  the  web  thickness  is  removed  due  to  yielding  in  shear  and 
with  the  use  of  von  Mises  yield  criterion,  the  remaining  web  thickness  used 
to  carry  normal  stresses  is 


WT  "  'w/1  "  3(s4V)2  <12> 

»       T  w  y 

When  V  is  equal  to  zero  the  special  case  of  pure  bending  occurs.   In  this 

case,  the  secondary  moment  due  to  shear,  M__,  is  equal  to  zero  and  w„ 

equals  t  . 
w 

Because  no  shear  is  applied  to  the  bottom  steel  tee,  it  provides  a 
constant  axial  tensile  force,  P  ,  throughout  the  low  shear  case  (Fig.  6) 

PB  =  (twSB  +  bt)Fy  <13> 


Force  P  has  a  corresponding  compressive  force  in  the  concrete  slab.  The 
thickness  of  the  concrete  slab  required  for  P  is  assigned  starting  from 
the  top  of  the  slab  and  is  determined  by 

PB 

(14) 


PB   b  F 
c  c 

The  forces  in  the  bottom  tee  and  concrete  slab  combine  to  give  the  primary 

moment.   To  find  this  moment,  the  distance  between  the  centroids  of  the 

two  forces  must  be  found.   From  Fig.  6,  the  distance  from  the  top  edge  of 

the  opening  to  the  line  of  action  of  the  force  in  the  concrete  slab  is 

yc  =  sT  +  t  +  c  -  %cpB  (15) 

while  the  distance  from  the  bottom  edge  of  the  opening  to  the  line  of 
action  of  the  force  in  the  bottom  tee  is 

%twSB2  +  bt(sB  +  kt) 

yB  =  1  sR  +  bt <16> 

w  B 

The  lever  arm  of  these  forces  is 

dc  =  yc  +  2h  +  yB  (17) 

thus  the  primary  moment  is  defined  as  the  product,  P  d  . 

Be 

There  are  two  cases  to  consider  in  the  low  shear  analysis  of  the  top 
steel  tee  -  concrete  slab  section  shown  in  Fig.  7  after  the  portion  of 
the  slab  due  to  the  primary  moment  is  removed.   These  are  Case  IA  in 
which  all  the  remaining  slab  in  Fig.  7  is  used  and  Case  IB  in  which  only 
part  of  the  slab  is  used.   The  location  of  the  NA-,  in  the  flange  or  the 
slab  of  the  section  in  Fig.  7  determines  at  the  outset  which  case  applies. 
To  determine  this  location,  the  axial  force  capacities  of  the  slab  with 
thickness 


cr  =  c  -  c 


PB  (18) 
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and  the  steel  tee  are  required.   They  are  respectively,  (Fig.  7) 

P    =  b  c  F  (19) 

ycr    c  r  c 

and 

PyT  =  (STWT  +  bt)Fy  (20) 

If  P  is  less  than  P  T,  then  the  NA_  is  in  the  flange.  Referring 
to  Fig,  8,  the  distance  to  the  NAp  in  the  flange  is  found  by  setting  the 
forces  above  and  below  equal  to  each  other  resulting  in 

STWT       bccrFc 

y m  st  +  *  -  -ir +  isP  (21) 

y 

Now  the  total  moment  capacity  of  this  section  by  summing  the  moments  about 
the  NAp  is 

Mcap  =  Vr(sT  +  t  -  y  +  %cr>Fc  +  ^V?  "  %ST} 

+  %b(y  -  sT)2  +  %b(sT  +  t  -  y)2]Fy  (22) 

When  a  non-zero  shear  is  imposed,  a  certain  portion  of  the  top  steel 
tee  is  assigned  a  moment  due  to  shear 

MvT  "  V  (23) 

This  shear  moment  is  assigned  to  the  extreme  top  and  bottom  edges  of 
the  steel  tee  moving  inward  and  is  restricted  by  the  location  of  the  NAp 
shown  in  Fig.  9.   The  portion  of  the  flange  above  the  NAp  is 

ty  =  sT  +  t  -  y  (24) 

and  a  depth  of  web 

btV 


SV  ~   w 
T 


(25) 


is  found  such  that  the  area  of  the  flange  above  the  NA^  is  equal  to  the 
area  of  the  web  corresponding  to  the  depth  s  .   If  s  is  less  than  s  as 


11 


shown  in  Fig.  9a,  then  the  distance  between  the  centroids  of  the  two 
forces  is  s _  +  t  -  %tv  -  %s„,  and  the  maximum  M^T  allowed  is  the  force 
times  its  lever  arm 

^max  =  bVsT  +  t  '  %tV  "  %sV)Fy  (26) 

When  sv  is  greater  than  s„  (Fig.  9b),  the  bottom  portion  of  MVT  goes 

into  the  flange  a  thickness 

-s  w  +  bt^ 
H.   "      b  (27) 

Summing  moments  about  the  NA_  gives 

"vmax  "  fSTWT^  "  ^  +  W   +  btVw^  '  *Vw  "  V]Fy    (28) 

In  both  cases  (s  greater  than  or  less  than  s„) ,  if  M^  is  less  than 

M„   ,  then  the  moment  due  to  eccentricity  is 
vmax'  ' 

Me  -  Mcap  -  "vi  <29> 

and  the  total  moment  capacity  of  the  beam  with  the  web  opening  is 

M  =  PBdc  +  Me  (30) 

When  VL      is  greater  than  M^   ,  part  of  the  slab  is  "thrown  away"  and 
Case  IB  is  encountered. 

Case  IB  with  the  NAp  in  the  slab  also  occurs  when  P    is  greater 
than  P   (Fig.  7).   This  second  major  breakdown  of  the  low  shear  case 
has  two  further  divisions  -  if  s   (as  described  previously)  is  less  than 
or  greater  than  s  . 

When  s  is  less  than  s„  as  in  Fig.  10a,  knowing  that  the  areas  in 
the  web  and  flange  must  be  equal,  the  thickness  of  the  flange  used  for 
M^  is 

SVWT 
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Using  the  force  and  lever  arm,  >L-T  becomes 

Mvt  =  SVWT(ST  +  t  ~  %tV  '  %SV)Fy  (32) 

but  is  also  equal  to  V  a.   Setting  these  two  equations  equal  and  substi- 
tuting for  tv  gives 

C%+5)SV2-  (8T  +  t)8V  +  ^jT  "  °  (33) 

T  y 

This  quadratic  equation  can  be  solved  for  s  ,  after  which  ty  can  be 
determined  from  Eq.  31.   Now  the  remaining  portions  of  the  web 

sp  =  sT  -  sv  (34) 

and  the  flange 

tp  =  t  -  tv  (35) 

are  used  to  find  the  axial  tensile  force  component  of  M  which  is 

PT  =    CspwT  +  btp)Fy  (36) 

An  equal  force  is  assigned  in  the  slab  starting  down  at  the  point  where 
cD  stops  until  the  thickness  as  given  by 

ID 

PT 

(3-7) 


PT   b  F 
c  c 

is  reached.   Summing  the  moments  of  these  two  forces  about  the  HA^    (which 

is  at  the  bottom  of  the  slab  being  used)  gives 

Me  =  %bccPT2Fc  +  [sPWT(cr  "  CPT  +  C  +  *P> 

+btp(cr  -   cpT+  tv+%tp)]Fy  (38) 

When  s  is  greater  than  sT,  the  bottom  portion  of  M^T  goes  into 
the  bottom  of  the  flange  as  in  Fig.  10b.   The  thickness  of  flange  above 
line  XX  on  the  top  tee  steel  section  now  becomes  by  setting  the  forces 
above  and  below  line  XX  equal 
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s  v 

*v  =  -T  +  Hw  (39> 

Summing  moments  about  the  line  XX  gives 

M^  =  [sTwT(t  -  tj  +  %sT)  +  %btv2  +  btVw(t  -  tv  -  %tVw)]Fy  (40) 

Equating  Eqs.  23  and  40  and  substituting  for  t^  results  in 

2 
btVw2  +  (STWT  "  bt)tVw  '  STWT(t  +  V  +  "f  +  "^"  =  °   (A1) 

y 

which  can  be  solved  for  t   .   Knowing  t..  ,  tv  is  found  by  Eq.  39  and  the 
thickness  of  the  flange  assigned  for  the  axial  force,  P  ,  is 

h  =  c  -  h  ~  ^w  (42) 

The  magnitude  of  the  axial  force  is 

PT  =  btpFy  (43) 

and  the  corresponding  force  equal  to  it  in  the  slab  has  thickness  cpT 
as  determined  by  Eq.  37.   The  moment  due  to  eccentricity  is  found  by  sum- 
ming the  moments  about  the  NA^  which  gives 

Me  =  %cpTPT  +  (cr  -  cpT  +  tv  +  %tp)PT  (44) 

In  both  cases  when  the  slab  is  not  completely  used,  the  total  plastic  moment 

capacity  is  given  by  Eq.  30. 

High  Shear  Solution.   The  second  major  case,  Case  II,  is  called  high 

shear,  in  which  part  of  the  total  shear  goes  to  the  bottom  tee  and  all 

the  top  tee  capacity  is  utilized  to  resist  V  and  M^T-   Because  the 

capacity  of  the  top  tee  is  used  entirely  for  V_  and  M™,  M  is  zero 

T      vT   e 

throughout  Case  II.   To  find  the  capacity  for  V  and  M^  of  the  top  tee, 
a  trial  and  error  method  is  applied  using  four  equations.  The  first  is 
the  expression  for  w  as  given  by  Eq.  12.   The  second  equation,  referring 
to  Fig.  11,  gives  the  thickness  of  the  flange  below  the  NA_  of  the  top 
steel  tee  as 
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-s  w  +  bt 
Sc  =  -^ (A5) 

Equation  23  is  the  third  equation  required,  and  the  last  one  is  found  by 

summing  moments  about  the  NA^  in  Fig.  11 

"VTI  =  [sTWT(tx  +  %ST}  +  %btx2  +  %b(t  "  V^y  (46) 

Assuming  a  value  of  V  ,  M^  and  M^i  are  calculated  and  compared  and  V 
is  adjusted   until  they  are  equal,  giving  the  capacity  of  the  top  tee 
for  V  and  M^T.   These  values  of  V„  and  M^T  are  constant  throughout  the 
high  shear  case.   With  the  shear  assigned  to  the  top  tee  known,  the 
shear  assigned  to  the  bottom  tee  is 

VB  =  V  -  VT  (47) 

and  the  moment  due  to  shear  in  the  bottom  tee  is 

^B  =  V  <48> 

Because  the  bottom  tee  now  has  shear  assigned  to  it,  it  has  a  reduced  web 
thickness 


WB  =  tvf-  3fcTF-)2  (49) 

t  B  w  y 

At  this  point,  the  treatment  of  the  bottom  tee  is  very  similar  to 

that  of  the  top  tee  in  the  low  shear  case  where  the  NA_  of  the  top  tee  - 

remaining  concrete  slab  section  was  in  the  slab.   The  calculations  are 

the  same  for  the  bottom  tee  as  the  top  tee  in  both  cases  (s  greater  than 

or  less  than  s  )  to  the  point  where  the  portions  of  the  tee  used  for 

the  axial  force  P  are  found. 

D 

When  s  is  less  than  sT,  the  axial  force  is  (Fig.  12a) 

PB  =  (spwB  +  btp)Fy  (50) 

The  corresponding  axial  force  in  the  concrete  is  assigned  to  the  slab 
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starting  at  the  top  and  having  thickness  c-,,  as  given  by  Eq.  14.  The 

CD 

distance,  y  ,  from  the  top  edge  of  the  opening  to  the  line  of  action  of 

the  force  P  in  the  concrete  is  expressed  by  Eq.  15  and  the  distance  from 

the  bottom  edge  of  the  opening  to  the  centroid  of  the  force  P  in  the 

B 

bottom  steel  tee  is 

%sp2wB  +  btpCsp  +  %tp) 

?b  =  — SpwB ;  btp +  sv  '     <51> 

The  moment  arm,  d  ,  of  the  forces  is  determined  by  Eq.  17,  and  is  used  to 
find  the  total  plastic  moment,  which  is 

M  -  PBdc  (52) 

because  M  is  zero, 
e 

In  the  other  case  of  s  being  greater  than  sT,  the  axial  force  is 


(Fig.  12b) 

PB  =  btPFy  <53> 

Again  the  same  force  in  the  concrete  is  assigned  starting  at  the  top  of 
the  slab  and  having  thickness  c  ,  which  is  calculated  from  Eq.  14.  The 
distance  y  to  the  line  of  action  of  the  force  P,.  in  the  concrete  from  the 

c  JJ 

top  edge  of  the  opening  is  given  by  Eq.  15,  while  the  distance  from  the 

bottom  edge  of  the  opening  to  the  centroid  of  the  force  P_  in  the  bottom 

B 

steel  tee  is 

yB  =  SB  +  *V*  +   %tP  (54) 

The  moment  arm  d  of  the  two  forces  is  determined  by  Eq.  17,  and  the  total 
moment  capacity  as  before  is  found  using  Eq.  52. 

Calculation  of  Interaction  Diagrams 

This  section  presents  the  sequence  of  calculations  used  in  developing 
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a  shear-moment  interaction  diagram.  A  broad  view  of  the  entire  sequence 
with  all  cases  will  be  presented  first,  with  the  details  of  each  individ- 
ual case  considered  later. 

Figure  13  is  the  overall  flow  diagram  of  the  procedure  followed  in 
developing  an  interaction  diagram.   First,  after  input  data  is  read, 
reference  values  for  a  composite  beam  with  known  dimensions  and  material 
properties  are  calculated.   One  limit  set  on  the  solution  at  the  outset 
is  that  the  total  axial  force  capacity  of  the  bottom  tee,  P  ,  must  be 

less  than  the  total  axial  force  capacity  of  the  concrete  slab,  P   .  This 

yc 

limit  is  used  since  a  composite  beam  with  the  force  P^  greater  than  the 

a 

force  P   is  an  impractical  case,  and  therefore  not  considered  here. 

If  P  is  less  than  P   ,  the  input  and  reference  values  are  printed, 
after  which  the  total  shear,  V,  (V  =  V  in  Case  I)  is  initialized  to  zero. 
The  value  by  which  the  total  shear  is  incremented  is  1.0  and  is  labeled 
V.   .   Later,  as  the  interaction  diagram  is  developed,  its  slope  becomes 
steeper,  requiring  a  smaller  increment  of  shear,  i.e.,  V.   =  0.1. 

At  this  point  a  program  control,  "check",  is  also  set  equal  to  zero. 
When  "check"  is  equal  to  zero,  a  further  decision  is  needed  before  going 
to  Case  IA  or  IB.  When  Case  IB  is  used  once,  "check"  is  set  equal  to  one, 
so  that  the  solution  process  returns  to  Case  IB. 

The  next  decision  deals  with  the  total  axial  force  capacities  of  the 
top  steel  tee  and  the  remaining  concrete  slab  (thickness  c  ) ,  which  are 
P  „  and  P   ,  respectively.   Details  of  this  decision  step  were  discussed 
in  the  previous  section.   After  this  decision,  the  solution  continues  to 
either  Case  IA  or  Case  IB,  both  of  which  are  shown  in  more  detail  in 
Figs.  IA  and  15,  respectively. 

At  the  end  of  either  case,  the  required  output  for  the  interaction 

diagram  is  printed.   The  value  of  shear  is  incremented  by  V.   and  the 

J      mc 
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new  shear,  V,  is  compared  with  the  total  allowable  shear  on  the  top  web 
stub,  V  .   If  the  value  of  shear  is  less  than  V   ,  then  the  process  is 
repeated  in  the  appropriate  case  giving  more  coordinates  for  the  inter- 
action diagram.   The  solution  is  stopped  if  V  is  greater  than  V   ,  since 
it  is  not  applicable  to  failure  in  shear. 

Case  IA  or  Case  IB  will  eventually  give  way  to  Case  II.  Figure  16 
is  a  detailed  flow  chart  of  the  solution  process  within  Case  II.  At 
the  end  of  Case  II,  data  for  the  interaction  diagram  is  printed  after 

which  the  shear  is  increased  by  V.   ,  which  is  now  0.1.   The  value  of  the 

mc 

shear  on  the  bottom  tee,  V_,  is  now  found  and  compared  with  the  total 

B 

shear  the  bottom  tee  stub  will  allow,  V   .   If  the  shear  force  V_  is  less 

yB  B 

than  V  ,  then  Case  II  is  repeated.   If  V  is  greater  than  V  ,  this 

solution  is  not  applicable  and  the  calculations  cease.  At  the  end,  enough 

coordinates  will  have  been  computed  to  plot  the  entire  interaction  diagram. 

Figure  14  shows  the  steps  involved  within  Case  IA,  all  of  which  have 

been  discussed  earlier  except  for  the  decision  of  whether  M  is  greater 

than  zero.  M  must  be  greater  than  zero  in  Case  IA  by  definition,  and 

if  it  is  not  Case  II  takes  over.  At  the  end  of  each  cycle  through  Case 

IA,  the  coordinates  of  the  interaction  diagram  are  computed. 

Case  IB  (Fig.  15)  is  activated  when  P  „,  is  less  than  P    or  M.T 

yT  ycr     Tmax 

is  less  than  M^.   The  value  of  "check"  is  changed  to  equal  1.0  so  that 
the  Case  IA  is  by-passed  through  the  remainder  of  the  solution.   The 
terms  AgV,  BgV,  CgV,  and  QgV  deal  with  the  quadratic  equation  for  s 
CEq.  33).  AgV,  BgV,  and  CgV  are  the  coefficients,  and  Q   is  the  portion 
under  the  square  root  of  the  quadratic.   If  Q   is  less  than  zero,  an 


imaginary  number  results,  so  the  solution  is  directed  to  solve  for  1 

V' 

in  a  manner  similar  to  that  for  s  .   If  Q    results  in  an  imaginary 

number,  the  solution  is  switched  to  Case  II.   If  either  s„  or  t   are 

V     Vw 


w 
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found,  the  remaining  calculations  are  performed,  and  coordinates  for  the 
interaction  diagram  are  computed.  Again,  a  check  for  M  is  made  in 
Case  IB  similar  to  that  in  Case  IA. 

Case  II  (Fig.  16)  occurs  when  M  is  less  than  or  equal  to  zero, 
or  when  Q  „  is  less  than  zero.  At  the  beginning  M  is  set  equal  to 
zero,  the  bottom  shear  to  top  shear  ratio  is  set  equal  to  zero  and  the 
value  of  shear  increment,  V.   ,  is  changed  to  0.1  for  reasons  given 
earlier.  With  the  given  shear  ratio,  V_  and  V„  are  found  and  the  moments 
My„   and  My„1  are  computed  and  compared.   Adjustments  are  made  to  the 
shear  ratio  until  M__  and  M--T1  are  equal.   Then,  as  in  Case  IB,  calcu- 
lations and  decisions  are  made  concerning  Q  TT  and  CvTT  .   If  Q  XT  is  less 

sv     xtVw      ^tVw 

than  zero,  the  solution  terminates.  Again  calculations  are  made  if 

values  for  s„  or  t„  are  found,  and  the  last  of  the  coordinates  for  the 
V     Vw 

interaction  diagram  are  determined. 
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TYPICAL  RESULTS  AND  DISCUSSION 

Interaction  Diagrams 

The  computer  solution  which  is  shown  in  Appendix  III  follows  the 
flow  diagrams  discussed  in  the  previous  chapter,  and  results  in  a  shear- 
moment  interaction  diagram  as  in  Fig.  17.   This  diagram  is  the  predicted 
failure  envelope  for  a  specific  beam  of  known  dimensions  and  material 
properties.   Shear  and  moment  are  non-dimensionalized  by  the  total  shear 
capacity  of  the  gross  web  section,  Vp,  and  the  total  plastic  moment 
capacity  of  the  gross  section,  M_,  ,  respectively.   For  any  given  set  of 
loading  conditions  and  opening  location,  the  theoretical  failure  load 
can  be  determined. 

As  indicated  in  Fig.  17,  two  possibilities  for  the  top  portion  of 
the  curve  were  investigated  based  on  two  different  methods  of  distribut- 
ing the  moment  due  to  shear  in  the  top  tee.   For  the  bottom  curve, 
Distribution  I,  the  moment  due  to  shear  was  assigned  at  the  top  of  the 
tee  section  as  shown  in  Fig.  18a.   The  interaction  diagram  from  this 
distribution  had  a  rather  sharp  downward  curve  at  the  beginning.   For 
Distribution  II  (top  curve)  the  moment  due  to  shear  was  assigned  at 
opposite  ends  of  the  top  steel  tee  (Fig.  18b) ,  resulting  in  a  higher 
moment  capacity  initially,  but  ending  with  a  slope  discontinuity  as  the 
two  curves  meet  at  the  end  of  Case  I.   Because  Distribution  II  gives 
a  higher  moment  capacity,  and  it  is  consistent  with  the  distribution 
assumed  in  the  bottom  tee,  it  was  adopted  for  this  analysis. 

The  slope  discontinuity  in  the  interaction  diagram  appears  to  be 
related  to  the  assignment  of  the  moment  due  to  shear  in  both  steel  tees. 
In  Case  I  the  total  moment  capacity  is  composed  of  the  primary  moment, 
which  is  constant,  and  the  moment  due  to  eccentricity,  M  ,  which  varies. 
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Because  the  primary  moment  is  constant  it  will  not  bring  about  a  change 

in  the  rate  of  decrease  of  the  total  moment  in  the  interaction  diagram, 

whereas  M  will.   The  change  in  M  is  brought  about  by  several  factors, 

the  first  of  which  deals  with  web  thickness.  As  shear  is  added  in  equal 

increments,  the  change  in  web  thickness  should  be  at  a  constant  rate 

thus  giving  a  constant  rate  of  change  in  the  interaction  diagram.  A 

second  factor  is  the  change  in  the  moment  arm  of  M  .  At  the  concrete 

end,  the  arm  would  be  increasing  as  less  concrete  is  used  for  larger 

shear  loads,  while  the  end  in  the  steel  will  become  shorter.   The 

concrete  is  not  "thrown  away"  faster  than  the  centroid  in  the  steel 

moves,  so  the  moment  arm  for  M  decreases  at  a  slight  rate  as  shear  is 

increased.   Since  the  magnitude  of  M  gets  smaller  as  its  moment  arm 

gets  smaller,  no  considerable  change  would  occur  in  the  slope  of  the 

interaction  diagram.   The  final  factor  deals  with  the  rate  at  which  area 

of  steel  is  used  for  M^  (or  Mr,B)  as  shear  is  added.  At  first,  a  small 

portion  of  the  top  tee  is  required  for  M™  because  of  a  large  moment 

arm,  but  as  more  shear  is  added,  more  area  of  steel  is  used  in  each 

increment  because  of  decreasing  moment  arm  length  (Fig.  19).  This 

would  cause  M  as  well  as  the  total  moment  to  become  smaller  at  an 
e 

increasing  rate,  giving  an  increased  rate  of  change  in  the  slope  of 
the  interaction  diagram.   The  slope  reaches  its  steepest  point  at  the 
end  of  Case  I,  after  which  in  Case  II  the  bottom  tee  is  assigned  M„ 
in  the  same  manner  as  the  top  tee,  so  the  slope  is  fairly  flat  at  first 
but  later  gets  very  steep. 

Figure  20  shows  a  comparison  of  the  interaction  diagrams  for  a 
non-composite  beam  and  a  composite  beam.   Both  curves  are  for  the  same 
W  shape  and  have  the  same  material  properties  and  opening  dimensions. 
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The  plot  for  the  non-composite  beam  was  produced  using  a  computer 
program  developed  by  Scritchf ield  (9) .   Because  the  beams  have  unequal 
total  plastic  moment  capacities,  the  M/l'L  coordinates  for  the  non- 
composite  beam  have  been  multiplied  by  Mp/Mp  to  permit  a  comparison. 
Since  the  composite  beam  has  a  higher  M/Kp  value,  it  would  appear  to 
be  the  more  effective  section.  At  the  lower  end  of  the  interaction 
diagram  the  two  curves  coincide,  which  should  be  expected  since  it  was 
assumed  that  the  concrete  does  not  carry  any  of  the  shear  force. 

Effects  of  Varying  Key  Parameters 

A  series  of  interaction  diagrams  have  been  prepared  to  investigate 

the  effect  of  some  of  the  key  parameters.   In  this  parametric  study,  a 

W  18x50  beam,  F  =  36  ksi. ,  f  =  3.5  ksi.  and  a  slab  width  of  48  in. 
y  c 

were  adopted,  while  slab  thickness  and  opening  length,  height  and 
eccentricity  were  varied  one  at  a  time.   In  the  following  discussion, 
an  interaction  diagram  for  c  =  A  in.,  h  =  4.5  in.,  a  =  6.75  in.  and 
e  =  0  is  common  to  all  of  the  figures. 

When  the  slab  thickness  is  varied,  not  much  change  is  effected  in 
the  interaction  diagram  as  can  be  seen  in  Fig.  21.   For  each  larger 
thickness ,  the  moment  capacity  for  any  value  of  shear  force  is  increased 
because  of  longer  moment  arms  for  both  M  and  the  primary  moment,  but 
the  total  moment  capacity,  VL    ,  is  also  increased,  resulting  in  little 
variation  in  the  K/VL      ratio.   Because  Kp  does  not  increase  faster  than 
the  moment  capacity  as  larger  thicknesses  are  used,  the  smaller  thicknesses 
have  larger  M/M^  values.  All  curves  meet  at  the  same  value  of  shear, 
showing  that  the  shear  load  is  independent  of  the  slab  thickness,  since 
it  is  assumed  that  the  slab  carries  no  shear. 
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Figure  22  shows  the  variation  in  the  interaction  diagram  for 
changes  in  opening  length.  With  a  shear  force  of  zero,  all  the  curves 
have  the  same  M/M^  ratio,  which  shows  that  change  in  opening  length 
does  not  affect  the  moment  capacity  in  pure  bending.   The  longer  the 
opening  length,  the  less  shear  load  the  beam  will  withstand.   This 
occurs  due  to  the  fact  that  moments  due  to  shear,  M„  =  V  a  and 
My  =  V_a,  increase  with  opening  length,  thus  with  a  longer  opening 
the  steel  section  is  spent  more  quickly  as  shear  force  is  increased. 

The  effect  of  varying  opening  height  is  illustrated  by  the  inter- 
action diagrams  in  Fig.  23.   The  smaller  the  opening  height,  the  greater 
the  M/M^  ratio  will  be,  because  less  of  the  beam  cross  section  is  lost 
to  the  opening.   Similarly,  with  the  smaller  opening  height,  a  larger 
shear  force  can  be  applied  to  the  beam  since  more  of  the  cross  section 
is  left  at  the  opening. 

Figure  24  shows  the  effects  on  the  interaction  diagrams  due  to 
variation  of  opening  eccentricity  (positive  eccentricity  is  upward  and 
negative  eccentricity  is  downward) .  The  largest  positive  eccentricity 
gives  the  highest  initial  M/>L,  ratio.   This  ratio  is  high  because 
steel  that  is  in  the  bottom  tee  will  have  a  larger  moment  arm  than  if 
it  were  in  the  top  tee.   As  the  eccentricity  decreases,  the  solution 
remains  in  Case  I  longer  since  more  steel  is  available  in  the  top  tee 
to  resist  shear.   Curves  with  equal  but  opposite  eccentricity,  closely 
converge  toward  the  bottom  portion,  suggesting  that  the  shear  capacity 
of  the  beam  is  not  significantly  affected  by  the  direction  of  eccentricity. 

Comparison  with  Experimental  Results 

Two  tests  of  composite  beams  with  web  openings  have  been  performed 
by  Granade  (5).  An  interaction  diagram  for  the  beams  is  shown  in  Fig.  25, 
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and  the  experimental  ultimate  loads  are  also  plotted.   A  large  dis- 
crepancy exists  between  the  theoretical  and  experimental  values  of  the 
failure  loads.   There  are  several  factors  which  might  contribute  to 
this  discrepancy;  however  their  effects  are  uncertain  because  the 
test  conditions  are  not  described  fully. 

A  small  factor  to  consider  would  be  the  manner  in  which  the 
material  properties  of  the  steel  and  concrete  were  determined.   This 
factor  would  cause  only  minor  changes  in  the  interaction  diagram. 
Another  small  change  might  occur  from  the  method  of  loading  the  beam. 
If  a  dynamic  loading  process  were  used,  a  higher  ultimate  load  would 
occur  giving  a  higher  test  point  on  the  interaction  diagram.  A  static 
loading  process  would  give  a  lower  ultimate  load.  The  effect  of  strain 
hardening  on  the  test  results  could  have  a  significant  effect.   Since 
the  ultimate  strength  analysis  does  not  take  into  account  the  effects 
of  strain  hardening,  the  experimental  ultimate  loads  would  have  to  be 
adjusted  (3)  to  give  a  good  comparison  between  theory  and  experiment. 

A  final  factor  concerns  one  of  the  key  assumptions  made  in  the 
analysis  presented  in  this  report.  The  assumption  states  that  no  shear 
force  will  be  assigned  to  the  concrete  slab.   If  part  of  the  shear  force 
were  assigned  to  the  slab,  ultimate  loads  predicted  from  the  interaction 
diagram  would  be  much  higher. 
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CONCLUSIONS 

An  ultimate  strength  analysis  of  composite  beams  with  web  openings 
has  been  developed  based  on  McCormick's  method.   This  analysis  was  used 
to  make  a  comparison  with  a  non-composite  beam,  and  the  composite  beam 
was  found  to  be  more  effective.  Ultimate  loads  based  on  this  solution 
were  also  compared  with  those  observed  in  two  laboratory  tests.   The 
theoretical  results  were  found  to  be  very  conservative  in  their  pre- 
dictions of  the  strength  of  the  test  beams. 

The  effect  of  variation  of  certain  parameters  of  a  composite  beam 
were  studied  using  the  analysis.   Observations  from  this  study  are  as 
follows : 

1.  Changes  in  the  slab  thickness  do  not  affect  the  interaction 
diagram  to  a  large  extent. 

2.  The  longer  the  opening  is,  the  smaller  the  failure  load. 

3.  As  the  opening  is  made  deeper,  the  moment  and  shear  capacity 
decrease. 

4.  An  opening  with  the  highest  positive  eccentricity  has  the 
highest  moment  capacity. 
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RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

Further  study  is  needed  in  regard  to  the  slope  discontinuity  in  the 
interaction  diagram.   This  study  should  be  directed  toward  determining  if 
an  assignment  of  forces  can  be  made  such  that  the  slope  discontinuity  is 
removed.  Also,  the  assignment  of  shear  force  to  the  concrete  slab  should 
be  considered  in  future  analytical  work.   The  analysis  presented  in  this 
report  could  be  expanded  so  that  it  could  be  applied  to  composite  beams 
with  reinforcement  at  the  web  opening. 

More  experimental  tests  on  composite  beams  with  web  openings  would 
be  helpful  for  comparison  with  theoretical  work. 
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APPENDIX  II  NOTATION 


a  -  one-half  length  of  opening 

b  -  width  of  steel  flange 

b  -  width  of  concrete  slab 
c 

c  -  thickness  of  concrete  slab 

c,,,,  -  thickness  of  concrete  used  to  equal  axial  force  P,, 

c„  -  thickness  of  concrete  used  to  equal  axial  force  P 
Ps  n  ys 

c   -  thickness  of  concrete  used  to  equal  axial  force  P„ 

c  -  thickness  of  concrete  left  after  thickness  cp  due  to  P_  is 
subtracted  from  original  thickness  c 

d  -  depth  of  steel  section 

d  -  moment  arm  between  axial  force  in  bottom  tee  and  corresponding 
force  in  slab 

d  -  moment  arm  between  axial  force  in  top  tee  and  corresponding 
force  in  slab 

e  -  eccentricity  of  opening 

F  -  .85  f' 
c       c 

i 

f  -  compressive  strength  of  concrete  cylinder 

F  -  yield  stress  of  steel 

h  -  one-half  opening  depth 

M  -  total  moment  capacity  of  beam  at  centerline  of  opening 

M    -  total  moment  capacity  of  top  tee-concrete  slab  (c  )  section 
cap  r 

M  -  moment  due  to  eccentricity 

e  ' 

VL   -  total  moment  capacity  of  non-composite  beam  without  opening 
Mpc  -  total  moment  capacity  of  composite  beam  without  opening 
Myg  -  moment  due  to  shear  in  bottom  tee 
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My    -  maximum  M^  allowed  in  top  tee  due  to  location  of  NA_ 

M^  -  moment  due  to  shear  in  top  tee 

My   -  value  of  My™  for  any  value  of  shear  by  £  Moments  -  used  to 
compare  with  value  My^. 

P  -  axial  force  in  bottom  tee  which  contributes  to  primary  moment 
B 

P„  -  axial  force  in  top  tee  which  contributes  to  M 
T  e 

P   -  total  axial  force  capacity  of  concrete  slab 

P    -  axial  force  of  concrete  slab  remaining  after  c_,_  removed 
ycr  6       PB 

P   -  total  axial  force  capacity  of  steel  section  at  opening 
P  _  -  axial  force  of  top  steel  tee  with  web  reduced  for  shear 
s  -  depth  of  web  section  in  bottom  tee  at  opening 

D 

s^  -  depth  of  web  assigned  to  axial  force  P  or  P 

r  B      T 

s  -  depth  of  web  section  in  top  tee  at  opening 

sv  -  depth  of  web  assigned  to  axial  force  component  of  M^  or  My 

t  -  steel  flange  thickness 

t  -  thickness  of  flange  assigned  to  axial  force  P  or  P 
r  B      T 

t  -  thickness  of  top  steel  flange  below  NAp  of  composite  beam 
without  opening 

t„  -  thickness  of  outside  edge  of  flange  assigned  to  M_  or  K 

tyw  -  thickness  of  flange  adjacent  to  web  assigned  to  M^  or  MVR 

t  -  steel  web  thickness 
w 

t  -  thickness  of  top  flange  in  tension  below  NA^ 

V  -  total  shear  applied  to  composite  beam  web  with  opening 
V  -  shear  assigned  to  bottom  tee 
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V  -  total  shear  capacity  of  web  of  steel  section  with  no  opening 

V  -  shear  assigned  to  top  tee 

V  -  total  shear  capacity  of  web  of  bottom  tee  section  at  opening 

V  -  total  shear  capacity  of  web  of  top  tee  section  at  opening 
w  -  reduced  web  thickness  for  bottom  tee 

w  -  reduced  web  thickness  for  top  tee 

y  -  distance  from  bottom  of  web  of  top  tee  section  to  the  NAp  of 
top  tee-concrete  slab  (c  )  section 

y„  -  distance  from  top  of  web  of  bottom  tee  to  centroid  of  portion 
assigned  to  axial  force  PR 

y  -  distance  from  bottom  of  web  of  top  tee  to  centroid  of  slab 
thickness  cpR  used  to  resist  force  PR 
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APPENDIX  III  COMPUTER  PROGRAM 


SJOB 
1                          REAL    M.MCJP.MF.Mf'PC.MPC.MVe.MVMAX.MVT.HVTGNE 
_2 PFADfS.I)    NBK 


3                     1    FORMATM5) 
*  DO    2000    J=l tNBK 

_5 eHFCK*CHFCKA=CHECKB*CHECKC«r.HFr.Kn  =  CKFCKF*CHECKF-CH£CK&-0 


6  P.EADI5.2I     E.O,T.TU,BC.E.H.A,FY,C,FPC 

7  ?   F0RMATI5F7.3.6F6.2) 

_B «;T.n/?.-'--n-T 


9  SB  =  D/2.  +  E-l--T 

10  VYT=ST*TW*FY/SQRTI3.) 

JJ VYR'SB»TH»FY/S0RT13.) 


12  VP=Tw*(D-7.*T»*FY/SCPT(3.) 

13  FC=0.E5*FPC 

JA PYC=FC»aC«C 


15  PYS«FY*(2.*P*T*(D-2.*T)*TW) 

16  IFIPYC    .LT.    PYS1G0    TO    3 

_LZ rPS*PYS/tFr«PC) 


18  MPC=FC*(BC»CPS**2/2.)  +  PYS*(D/2.*C-CPS) 

19  CO    TO    4 

_2fl 3   TT  =  lFr.»BC*r-FY»TW»(D-7.*TI)/(?.»FV«B) 


21                        MPC«FC*(BC*C*(C/2.+T-TT)  )*FY*(  I  T-TT  )**2*3/2.+ ITT)**2«B/2.+Tli*  f  C-2. 
C*T)*(D/2.-T*TT )+T*B*(D-3.*T/2.*TT>) 
JLZ 4    PB=FY»1B»T*SP»TH) 

23  IFIPd    .LT.    PYCIGO    TO    6 

24  HRITE(6.5I 

_25 S    Fr)R^T(7X.141HSINCF    THF    AXIAL    YIFLO    FORCE     IN    T)- F    BCTTCH    TEE    IS    GRE 

IATER    TI-AN   THE    AXIAL    FORCE    CAPACITY    CF    THE    CONCRETE    SLAB,    THIS    SOLU 
1TICN    IS    NOT    APPLICABLE.) 

_2A cn   Tfi   ?,vifi 

27  6    KRITFI6.7I     E.  D.T  ,TW.  BC.  E  .  H,  A.  FY  ,C  ,FPC 

28  7    FORMAT<lH0.<.X,lHE.9X,lhC.9X,lHT.9X.2HTW,7X.2HBC.9X,lHE.9X.lHH.9X,l 
CHA.BX.2HFY.CX.1HC.RX.3HFPC.//11F1C.3.//1 

29  WRITE(6,8>     ST , SB .VYT .VYB.VP.MPC 

30  8   FCRM4T(lH0.I0X.?HST,13X,2HSB.t2X,3HVYT,l2X,3HVYE,13X,2HVP,llX»3HHP 

r.r..//feFis.3.//  i 

31  WRITEI6.9I 

32  9   F0RMAT(lH0,10X,lHV.13X.5HVB/VT,llX,2HDC,l2X,lHf>  ,15X,2HME, 12X.4HV/V 
CP.11X.5HH/HPC) 

33  V=VT»0 

34  VBTO 


_35 VIKC1.0 


36  10   WT"TW*SORTI  l.-3.*(VT/(ST*TW*FYI  »*»2I 

37  CPB»r>3/!0C»FC> 

_3J VC*ST»T*C-CPB/2. 


39  YB=(l./(SB*TW*B*T))*(SB**2*TW/2.+B*T*(SB»T/2.)) 

40  DC*YR+H*2.*YC 

41 CR»C-CPB 


42  PYCR«FC*BC»CR 

43  PYT=FY*(B*T+ST*WT) 
44 HVT«VT»A 


45  IFlChECK    .EO.    1.)    GO   TO    27 

46  IFICHFCKA    .EC.    1.160   TO    15 
-42 IFIPYCR     .LT.    PYTIGC    TO    13 


4R  II    WRITEI6.12) 

49  12    FflRMAT(lH0,40HPYCR    IS    GREATER    THAN    PYT.    GO    TO    CASE    IB.) 

_5_2 GP  TQ   ^7 


51  13    WRITFI6.14) 

52  14    FORMAT! 1H0.37HPYCR    IS    LESS    THAN   PYT,    GO    TC    CASE    IA.) 
-51 CHECKA'l 
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54  15   Y=ST*T/2.-$T*WT/(7.*B>*(FC/FY)*tCR*BC/(2.»BI) 

55  MCAP=FY*IST*KT*(Y-ST/2.)*B/2.*(Y-ST)**2«-B/2.*(SH-T-Y)**2)»CR*BC*IS 

n  +  T-v«ca/?.i»FC 

56  TV=ST+T-Y 
5T  SV-TV*E/tfT 

_Sfl IF(SV    .  ST.    ST  ICC   TC    17 


59  MVMAX*TV*B*IST*T-TV/2.-SV/2.l»FY 

60  irjCHECKB    .EC.    l.)GC   TO    19 
_6J WR1TCI6.U) 


6?                 16    FtlRMATUH0.4X,24HME    FXTENDS    INTC    THE    WEB.) 
63                         CHECKB«1 
.64 GO    TO    19 


65  17    TVW«(TV**-ST*WT)/B 

66  MVKAX=(TV»*2/2.*P*TVW*B«-(T-TVW/2.-TV>*ST*WT*IY-ST/2.))*FY 

67 1FICHECKC  .EC.  1.1GC  TO  19 


66  WRITEI6.18I 

69 

70 


18  FCRMATI1H0,4X,2VHME  IS  CONFINED  TO  THE  FLANGE. I 
CHECKC-1 


71  19    IFIMVMAX    .GT.     MVT)GC    TO   21 

72  WRITEI6.20) 

73 20    rOPHAT(lH0.38HMVKAX     IS    LESS    THAN    MVT  .    GO    TO   CASE    IB.) 


74  GO    TO    27 

75  21    ME.MCAP-MVT 

76 22    IF  (WE    .GT.    0.)    GC    TO    24 


77  WRITEI6.23) 

78  23   FORMATUH0.36HME    IS    LESS    THAN    2ERC.    GO    TO   CASE    II.) 
79 GH    TO    33 


80  24   M=PB»DC*ME 

81  VVP»V/VP 
82 MMPC«M/»PC 


83  WP.ITEC6.25)    V . VBT ,DC . M.ME . VVP ,MMPC 

84  25   FCPMATI7F15.4) 

85 V»VT*V»VINC  


86  IFIVT     .LT.    VYT1G0    TO    10 

87  WRITEC6.26) 

88 26  FQRHAT(1H0.61HWHEN  VT  IS  GREATER  THAN  VYT ,  THIS  SOLUTION  IS  NOT  AP 

1PL1CABLE.) 
89  GC  TO  2000 
90 77  CHECK-1 

91  ASV*.5*WT/I2.*B) 

92  BSV«-(ST*T) 

93 CSV=V»A/IKT»FY) 

94  0SV«RSV**2-4.«ASV»CSV  '■ 

95  IFIOSV    .LT.    C.IGC    TO    29 

96 SV=l-eSV-SCPT  [QSV)  )/(2.*ASV)    


57  IF(SV    .GT.     STIGC    TC    29 

98  TV*SV*WT/B 

99 SP«ST-SV 


100  TP*T-TV 

101  PT«FY*(TP*0»SP»WTI 
10? CPT«PT/1BC«FC) 


103  MF  =  IBC*C<»T»»2/2.  )*FC*  t  SP*WT*  (CR-CPT«T»SP/2  .  )«-TP«8»  ICR-CPT»TV*TP/2. 

C)  )*FY 
104 1FCCHCCKD    .EC.    l.)CC    TC    22 


105  WRITEI6.28) 

106  2R    FORMAT  C1H0.4X.44HSV    IS    LESS    THAN    ST:    ME    EXTENDS    IMC    THE    WEB.) 
1C7 CHECKO*! ___ 

1 08  GO    TO    22  ~ 

109  29    ATVW«P 

110 BTVW»(ST»WT-B«T)  
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HI                        CTVW=-ST»UT*(T+ST/2.)*V*A/FY*!ST*t.T)**2/I2.*BI 
112                         0TVW=MBTVW»*?)-4.»ATVW»CTVW 
_LL2 IF10TVW    .GT.    O.ICO    TO    31 


114  WRITEI6.301 

115  30  FORMAT! 1H0 , 6OHS0UARE  ROOT  IN  CUADRATIC  FOR  TVW  IS  NEGATIVEt  GO  TO 
LEASE  U.) 


116  GO    TO    33 

117  31    TVH  =  !-RTVW-SORTICTVW) )/I2.*ATVWI 

1  18 TV*IST«HT1/P4-TVW 


119  TP«T-TV-TVk 

120  PT«TP*B*FY 
_L2J f.PT  =  PT/IBf»Fr> 


122  ME=CPT*PT/?.+!CR-CFT+TV*TP/2.l»PT 

123  IFCCHECKE    .EC.     l.)GO    TO    22 

.124 WRITE16.3?! 


125                32   FORMAT!  1H0.4X.52HSV    IS    GREATER    THAN    ST;    ME    IS    CCNFINEC    TO    ThE    FLAN 
1GE.I 
-L2A (-.HFCKF.I 


127          GO  TO  22 
12B       33  ME»0 
_1Z9 VBT.Q 


130  34   VT*V/I l.+VBTJ 

131  VB-V-VT 

_U2 iFtvT   .it.   vYTtr.r  Tr  35 


133  VBT=VBT*.301 

134  GG    TO   34 

135 3^   TFtvn    .i  t.   WBir.n  to  36 


136  VBT=VBT-.001 

137  GO    TO    34 

-L3J3 36   WT=TH*S0RT!l.-3.»fVT/IST«Th»FY)  I»»?l 


139  TX«(B*T-ST*WT 1/(2. *BI 

140  MVT=VT*A 

JLAJ MVT0NF*HT-TXl«»;>«B/?.+B«TX««2/2.»ST»HT»!TX+ST/7.  1>»FY 


142  IFIKVTCNE    .GT.    MVT)GC    TO    38 

143  VBT«VBT*.0001 
144 r.fl    TO    34 


145  37    VBT=VB/VT 

146  MVB*V8«A 

JU2 3P  wR»Tw»s:RTn.-i.*ivn/(«;R*Tw»FY)  >««?» 


148  ASV=.5+WR/(2.*8) 

149  ESV»-!SB«-T) 

151) CSV=VR«A/(WR»FVI 


151  0SV=BSV*»2-4.*ASV*CSV 

152  IFIOSV    .LT.    O.IGC    TO    40 
-132 SV*(-PSV-SCRT!CSH)/!?.»ASV> 


154  IFCSV    .GT.    SBIGO    TO    43 

155  TV«SV*WB/B 

-156 TP*T-TV 


157  SP=SB-SV 

158  PB«FY*!SP*WP+TP*B) 

159 YF=!TP«B«I  SP+TP/?.)+SP»*7*WB/2. ) / 1 SP»WB*TP« B) *SV 

160  IFICHECKF  .EO.  l.IGO  TO  44 

161  WRITEI6.39) 

-U>2 39    FORMAT!  1H0.4X.44HSV    IS    1FSS    THAN    SB;    PB    EXTENDS    I NTC    THE    WEB.) 

163  CHECKF=1 

164  GO  TO  44 

165 4T  ATVtt=B 


166  BTVW=(SB*WB-B*T) 

167  rTVH=-SRi-WB*IT*SR/2.  )»VB*A/FY+<  SB*Wfil**2/  !2.*Bl 
-1A8 OTtfy*  (  RTVU««7  I-4  .♦ATVW»rTtfU 


34 


169 

170 

171 

41 

IF  ICTVt.    .GT.    O.IGC    TC    42 
WRITE (6. 411 

FORMAT  I1HO.55HS0UARC    RCl.T     IM    TFE    CUAOKATIC    FOR 

TVV.    IS    NEGATIVE  .    ST 

17? 
173 

42 

OP.) 

r.c    TO    2  3  TO  • 

T\/w=(-BT\/Vi-SOr<TICTV.<)  l/I2.«ATVWJ 

174 
175 

176 

TV*<  Sa*W3)/Q*TVW 

TP=T-TV-TVW 

PR*TP*P*FY 

177 
178 
179 

Y3  =  Sft«-TVV.*TP/2. 

IF(CHEC<r.    .EO.    l.IGC    TC    44 

WRlTElt.431 

180 
181 

43 

FORMAT)  1HJ.4X.52HSV    IS    GREATER    THAN 
LGF.I 
CHECKG=1 

SB:   p a   is 

CCNFINEC    TC    THE    FLAN 

id? 

183 

194 

44 

VINt=  3.1 

cpp=pe/tec*FO 

YC=ST+T*C-CP8/2. 

1  35 
166 
}87 

DC=Y3»2.*h*YC 

M=PD*CC 

VVP=V/VP 

188 
189 
1 90 

45 

MNPC'P/.MPC 

WRITE (6. 451    V.VBT.CC.f rfE.VVP.MMPC 

F0RMATI7F15.4) 

191 
19? 

193 

V=V+VINC 

VS*V-VT 

IFIV9     .LT.    VYB)<",Q    TO    37 

194 
195 

46 

WRITE(fc.46) 

FCRMATUH0.61HHHEN    V8     IS    GREATER    THAN    VY3.    THIS    SCLUTICN     IS    NCT    AP 
1PLICA81E.)                                                                                                                                                   ■■■-    - 

196 

"]                  197 

198 

2000 

CCfaiNUE 
STOP 

FND 

iENTRY 

a 

7.500 

0                          T                          TW                    QC 
ia.300               0.570               0.358            43. COO 

E 

0.30C 

H                          A 

4.5CC               6.753 

FY 

36.030 

C                       FPC 
4.300              3.53C 

" 

SI 

3.930 

se                            VYT 
3.930                          29.243 

VY8 
29.243 

VP 
125.452 

MPC 
566C.773 

„ 

"                   pYCR 

V                                    VB/VT                                CC 
IS    CREATFP    ThAN    PYT,     CC    TC    CASE     IE. 

N 

*E 

V/\iP 

I'/fPC 

sv 

IS    ITSS    THAN    ST:    ME    EXTENDS    INTO    THE    WEB. 
0.0000                          0.0000                       20.4410 
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Fig.  1  Elevation  of  Composite  Beam  with  Web  Opening 
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Fig.  2   Section  of  Composite  Beam  with  Web  Opening 
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Fig.  3  Four  Hinge  Failure  Mechanism 


41 


<> 


T     T  - 


?T 

_PT 


s 


a.      Case   I     Pure   Bending    (V  =   0,   M  =   Pgd     +  Mg,   Mg  -  P^) 


B 
PT- 


(I 


fJV 


V 


s 
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Fig.  14  Flow  Diagram  for  Case  IA 
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Fig.  16  Flow  Diagram  for  Case  II 
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Fig.  18  Methods  of  Shear  Moment  Distribution 
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Fig,  20  Interaction  Diagrams  for  Composite  and  Non-composite  Beams 
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Fig.  21  Effect  of  Varying  Slab  Thickness 
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Fig.  22  Effect  of  Varying  Opening  Length 
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Fig.  24  Effect  of  Varying  Eccentricity 
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Fig.  25  Test  Results  from  Reference  5 
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ABSTRACT 

The  purpose  of  this  thesis  is  to  present  an  ultimate  strength 
analysis  of  composite  beams  with  web  openings.  With  the  use  of  this 
analysis  certain  variables  were  studied  and  the  following  conclusions 
were  drawn: 

1.  Changes  in  the  slab  thickness  do  not  affect  the  interaction 
diagram  to  a  large  extent. 

2.  The  longer  the  opening  is,  the  smaller  the  failure  load, 

3.  As  the  opening  is  made  deeper,  the  moment  and  shear  capacity 
decrease. 

4.  An  opening  with  the  highest  positive  eccentricity  has  the 
highest  moment  capacity. 

Theoretical  results  based  on  the  analysis  provide  a  very  conserva- 
tive prediction  of  the  strength  of  test  beams.   This  is  thought  to  be 
primarily  due  to  the  assumption  that  the  concrete  slab  does  not  carry 
any  shear  force. 


